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PROTON  DEPTH  DOSE  DISTRIBUTION:  3-D  CALCULATION 

OF 

DOSE  DISTRIBUTIONS  FROM  SOLAR  FLARE  IRRADIATION 


PURPOSE 


This  work  was  designed  to  estimate  the  relative  dose  distributions  to 
the  human  head  from  typical  solar  flare  irradiation  under  3  different 
circumstances: 

1.  Single  directional  radiation  incident  upon  a  fixed  head  (fixed 
field  irradiation). 

2.  Single  directional  radiation  incident  upon  a  head  rotating 
axially  (2-dimensional  (2-D)  rotation). 

3.  Omnidirectional  radiation  incident  upon  a  head  (3-dimensional 
(3"D)  rotation). 

This  work  differs  from  other  previously  reported  work  in  that  calculations 
are  made  relative  to  the  actual  3-D  shape  of  a  human  head,  rather  than 
estimates  based  on  a  cylindrical  shape.  The  calculated  isodose 
distributions  display  the  influence  of  the  head's  shape  for  each 
calculation  circumstance  just  listed. 


DATA  USED 


Three  sets  of  solar  flare  depth  dose  data  were  used  in  these 
calculations.  The  data  were  accumulated  from  measurement  of  1  solar  flare 
and  2  cyclotron  simulations  of  solar  flares.  The  3  sets  of  solar  flare 
data  are: 


1.  Simulated  solar  proton  event  of  10  July  1959.  (Approximately  1 0% 
dose  at  2-cm  depth)  Reference:  B.S.  Burton,  "Simulation  and  Measurement  of 
Solar-Flare  Proton  Exposures,"  June  1970,  SAM-TR-70-38. 


2.  Simulated  solar  proton  event  of  23  February  1956.  (Approximately 
1 0%  dose  at  5-cm  depth).  Reference:  B.S.  Burton,  "Simulation  and 
Measurement  of  Solar-Flare  Proton  Exposures,"  June  1970,  SAM-TR-70-38. 


3.  Depth-dose  spectrum  for  a  DMSP/F7  crossing  of  the  north  polar 
cap  at  the  peak  of  the  solar  particle  event  on  16  February  1 98-4  (10$  dose 
estimated  at  9.25-cm  depth).  Reference:  M.S.  Gussenhoven ,  E.G.  Mullen, 
R.C.  Filz,  D.H.  Brautigam,  and  F.A.  Hanser ,  "New  Low-Altitude  Dose 
Measurement,"  IEEE  Transactions  on  Nuclear  Science,  Vol.  NS-3^,  No. 3,  June 
1987,  pp.  676-683. 
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CALCULATION'  TECHNIQUE 


An  existing  clinical  radiation  therapy  treatment  planning  program  was 
used  to  calculate  the  single  fixed  field  isodose  distributions  and  the  2-D 
rotational  isodose  distributions.  A  3-D  dose  calculation  program  was 
written  to  calculate  the  3-D  isodose  distributions.  Key  points  of  the  dose 
calculation  program  for  all  3  oases  are: 

1.  All  isodose  distributions  are  displayed  relative  to  a 
normalization  dose  of  100  centigray  at  the  isocenter  in  the  absence  of  the 
head  phantom. 

2.  The  incident  beam  profile  is  assumed  to  be  uniform  (no  variation 
in  intensity  across  the  field  in  a  direction  perpendicular  to  the  incident 
ray) . 


3.  Effective  source  position  for  all  calculations  set  to  1,000  cm 
(this  simulates  a  planar  field  with  only  minimum  divergence). 

4.  The  head  phantom  is  assumed  to  be  of  unit  density  (no  correction 
for  increased  density  of  skull). 

5.  Tissue  absorption  effects  are  calculated  by  determining  the  ray 
path  length  through  the  head  phantom  to  the  calculation  point.  The  ray 
originates  at  the  effective  source  position  (1,000  cm  from  isocenter)  and 
continues  through  the  calculation  point.  The  typical  path  length  through 
the  head  phantom  to  isocenter  ranges  from  7  cm  to  12  cm,  while  the  path 
length  for  points  on  the  surface  of  the  phantom  is  zero  for  beams  incident 
on  that  surface  and  may  be  as  much  as  20  cm  to  that  same  point  for  beams 
incident  on  the  opposite  side  of  the  head. 

6.  An  inverse-square  correction  is  included  for  all  calculation 
points.  For  a  calculation  point  at  isocenter,  this  factor  is  1.0;  for  a 
point  10-cm  "upstream"  from  isocenter  relative  to  a  fixed  field,  this 
factor  is  (1,000/990)  squared,  or  1.01;  for  a  point  10cm  "downstream"  from 
isocenter  relative  to  a  fixed  field,  this  factor  is  (1,000/1,010)  squared, 
or  0.98. 


7.  Doses  for  the  single  fixed  field  and  2-D  rotations  are  calculated 
using  a  dose  grid  of  approximately  9,500  points,  while  doses  for  the  3-D 
rotations  are  calculated  using  a  dose  grid  of  approximately  11,900  points. 

8.  Angular  increment  for  the  2-D  rotation  is  set  to  3  degrees;  the 
angular  increment  for  the  3-D  rotation  is  set  to  20  degrees  and  the  entire 
rotation  is  simulated  by  14  arc  segments  equally  distributed  over  the  head. 
Total  dose  to  a  calculation  point  is  the  sum  of  dose  contributions  from  each 
increment.  The  individual  weight  applied  to  the  dose  from  each  increment  is 
(1/(total  number  of  arc  increments)). 

9.  Depth  dose  data  are  tabulated  in  2.5-mm  increments  to  a  depth  of 
10  cm.  For  depths  greater  than  10  cm,  calculations  are  based  on 
logarithmic  extrapolation  from  the  last  two  points  in  the  depth  dose  table. 
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VERIFICATION  OF  CALCULATIONS 


The  3~D  dose  calculation  program  was  intercompared  with  a  similar 
program  in  use  at  the  Joint  Center  for  Radiation  Therapy,  Harvard  Medical 
School,  Boston,  Mass.  This  intercomparison,  using  the  same  data  set  for 
the  head  phantom  contours  and  the  same  data  set  for  the  solar  flare  depth 
doses,  showed  that  the  isodose  distributions  resulting  from  calculations  on 
the  2  independently  developed  systems  agreed  to  within  the  thickness  of  the 
isodose  line.  Additionally,  the  depth  dose  distributions  corresponding  to 
2-D  rotation  were  calculated  for  a  10-cm  diameter  cylinder  and  compared 
with  the  rotated  depth  doses  reported  for  simulated  solar  proton  events  of 
10  July  1959  and  23  February  1956.  This  comparison  revealed  differences 
of  approximately  6$  in  the  calculated  surface  dose,  a  maximum  deviation  of 
1 0%  at  one  calculation  point  in  a  high  gradient  region,  and  deviations  of 
3%  or  less  for  depths  greater  than  1.5  cm.  We  referred  back  to  our  original 
data  interpolation  from  the  small  depth  dose  graphs  which  had  been  plotted 
on  semi-logarithmic  scale,  magnified  these  graphs  to  full-page  size  and 
repeated  the  data  interpolation  from  the  larger  graphs.  Based  on  this 
information,  we  located  interpolation  differences  of  up  to  approximately 
6%.  The  revised  depth  dose  data  set  was  entered  into  the  standard 
treatment  planning  program  data  base  and  the  calculations  were  repeated. 
Using  the  revised  data  set,  the  deviation  remained  6%  on  the  surface  (.501 
compared  to  .532),  with  a  maximum  difference  of  7.5%  (.151  compared  to 
.163)  at  a  depth  of  0.5  cm.  Between  depth  of  1.0  cm  and  isocenter  depth  of 
5  cm,  agreement  is  improved  to  within  2%.  Using  depth  dose  data  for  the 
Simulated  Proton  Event  10  July  1959,  plane-wave  calculations  of  surface 
dose  for  cylindrical  phantom  diameters  ranging  from  10  cm  to  20  cm  showed 
predicted  surface  dose  varying  from  52%  to  53%  of  the  incident  flux.  The 
2-D  rotation  calculation  on  the  phantom  skull  predicted  a  surface  dose  of 
M6 . 1  %  of  the  incident  flux.  This  calculation  suggests  that  the  isodose 
distributions  displayed  in  this  report  may  underestimate  the  maximum 
surface  dose  by  as  much  as  10^  (.^61  compared  to  .52),  while  relative  doses 
at  depths  below  the  phantom  surface  will  be  estimated  to  within  3%.  This 
result  may  reflect  a  potential  limitation  of  applying  the  isodose 
calculation  techniques,  which  interpolate  between  dose  grid  points,  to 
determination  of  dose  in  a  high  gradient  region.  For  the  Simulated  Proton 
Event  10  July  1959,  a  positional  shift  or  i  mm  introduces  a  \5%  change  in 
depth  dose.  Additional  work  is  needed  to  evaluate  this  difference  in 
surface  dose  predictions  with  calculation  technique. 


RESULTS 

Calculated  isodose  distributions  in  the  transverse  plane  intersecting 
isocenter  are  included  in  Figures  1-15  for  each  of  the  3  solar  flare  events 
for  fixed  field  irradiation,  2-D  rotation,  and  3“D  rotation.  For  the  2-D 
rotation, doses  are  also  reported  in  the  most  superior  plane,  corresponding 
to  the  crown  of  the  head.  Additionally,  isodose  distributions  are  recorded 
in  the  sagittal  plane  intersecting  isocenter  for  the  3”D  rotation.  In  all 
cases,  the  isodoses  are  reported  relative  to  the  dose  at  isocenter  (100 
centigray)  in  the  absence  of  the  head  phantom.  Wire-frame  models  of  the 
head  phantom  indicating  the  calculation  plane  and  the  rotation  plane  for 
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each  arc  segment  in  the  3“D  rotation  are  included  in  pages  22-36.  A 
comparison  of  dose  at  isocenter  vs.  dose  at  the  phantom  surface  is 
presented  for  each  solar  flare  event  in  Table  1 . 


EVALUATION  OP  RESULTS 

The  isodcse  calculations  for  the  3  solar  particle  events  reviewed 
show  that,  in  all  calculation  configurations,  the  maximum  predicted  dose 
will  be  on  the  exterior  surface  of  the  head  phantom.  The  relative  surface 
dose  to  the  phantom  increases  with  increasing  energy  of  the  solar  flare 
event;  this  is  attributed  to  the  greater  exit  dose  achieved  due  to 
increased  penetrability  of  the  solar  particles  with  the  increasing  energy. 
The  minimum  dose  within  the  phantom  incre-'es  with  the  energy  of  the  solar 
event.  For  a  given  solar  event,  the  minimum  dose  within  the  head  phantom 
is  approximately  the  same  for  both  2-D  and  3-D  rotation  calculations, 
indicating  that  the  average  depth  to  that  calculation  point  remains 
approximately  the  same  for  both  calculations.  This  depth  is  realistic  for 
the  head;  however  we  would  expect  to  see  a  greater  difference  between  2-D 
and  3“D  rotation  calculations  if  applied  to  the  thorax,  abdomen,  arms  or 
legs  due  to  the  elongated  shape  compared  to  the  head.  This  calculational 
technique  can  equally  well  be  applied  to  predictive  display  of  dose 
distributions  throughout  the  thorax  and  abdomen  and  can  be  used  to 
demonstrate  the  effect  of  tissue  inhomogeneities  such  as  lungs  and  bones  on 
the  resultant  doses  received  during  solar  events.  We  suggest  that  this  be 
investigated  in  a  follow-on  study  using  these  solar  flare  data  as  well  as 
other  representative  examples. 

These  calculations  suggest  that  the  superficially  located  organs  are 
at  greatest  risk  to  radiation  similar  to  the  3  cases  studied  here.  Thus, 
the  lens  of  the  eye  and  the  skin  would  receive  the  greatest  exposure. 
Because  of  the  greater  sensitivity  of  the  lens,  it  will  be  the  limiting 
organ  in  determining  necessary  shielding  and  exposure  time  limits. 
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TABLE  1.  SIMULATED  PROTON  EVENT  10  JULY  1959 
(CASE  1) 


Dose  at  Isocenter 

Dose  at  Surface 

Ratio 

Fixed  Field 

1 .0 

101.3 

2-D  Rotation 

0.8 

46.1 

0.017 

(superior  plane) 

3.55 

47.1 

0.075 

3-D  Rotation 

0.75 

27.8 

0.027 

SIMULATED  PROTON 

I  EVENT  23  FEBRUARY 

1956 

(CASE 

2) 

Dose  at  Isocenter 

Dcce  at  Surface 

Ratio 

Fixed  Field 

6. 

101 .11 

2-D  Rotation 

5.1 

50.3 

0.101 

(superior  plane) 

12.0 

54.0 

0.222 

3~D  Rotation 

4.9 

38.0 

0.128 

SIMULATED  PROTON  EVENT  16  FEBRUARY 
(CASE  3) 

Dose  at  Isocenter  Do3e  at  Surface 

1984 

Ratio 

Fixed  Field 

13. 

103.5 

2-D  Rotation 

11.1 

57.8 

0.192 

(superior  plane) 

23.9 

63.1 

0.379 

3-D  Rotation 

10.6 

55.6 

0.191 
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Figure  1.  Simulated  proton. event  10  July  1959.  Fixed  field  isodose  distribution 
(isocenter  transverse  plane) . 
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Figure  2.  Simulated  proton  event  10  July  1959.  Two-dimensional  rotation  isodose 
distribution  (isocenter  transverse  plane). 
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Figure  3.  Simulated  proton  event  10  July  1959.  Two-dimensional  rotation  isodose 
distribution  (superior  transverse  plane). 
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Figure  4.  Simulated  proton  event  10  July  1959.  Three-dimensional  rotation  iso¬ 
dose  distribution  (isocenter  transverse  plane) . 
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Figure  5.  Simulated  proton  event  10  July  1959.  Three-dimensional  rotation  iso¬ 
dose  distribution  (isocenter  sagittal  plane) . 
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Figure  7.  Simulated  proton  event  23  February  1956.  Two-dimensional  rotation 
isodose  distribution  (isocenter  transverse  plane). 
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Figure  11.  Solar  particle  event  16  February  1984.  Fixed  field  isodose  distri¬ 
bution  (isocenter  transverse  plane). 
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Figure  15.  Solar  particle  event  16  February  1984.  Three-dimensional  rotation 
isodose  distribution  (isocenter  sagittal  plane). 


The  following  wire-frame  graphics  displays  illustrate  the  3-D 
orientation  of  the  head  phantom  during  the  sequence  of  arc  rotation 
calculations. 


p.22  3-D  reconstruction  of  the  head  phantom  from  the  computerized 

tomography  scans.  The  top  two  views  demonstrate  the  phantom  contour 
outline  in  each  of  20  CT  slices.  The  bottom  two  views  show  the  wire¬ 
frame  reconstruction  of  the  head  phantom  1.3  viewed  from  the  side  and 
the  front. 


P-23 


p.24 

p.25 


p.26 

p.27 

p.28 

p.29 

p.30 

P-31 

P.32 

P-33 

P-34 

P-35 

P-36 


Head  positioned  -MO  degrees  relative  to  normal  incidence  of  radiation 
field.  Upper  left:  Lateral  view  of  head  as  rotated  -40  degrees. 
Vertical  line  is  plane  of  rotation  of  radiation  field.  Solid  contour 
line  is  intersection  of  transverse  calculation  plane  with  head 
phantom.  Lower  left:  Anterior  view  of  head.  Horizontal  line  is 
plane  of  rotation  of  radiation  field.  Line  at  40  degrees  below 
horizon  is  transverse  calculation  plane.  Lower  right:  Limits  of  arc 
•segment  for  this  orientation.  In  all  succeeding  displays  the  same 
general  description  applies.  Only  the  angle  of  the  head  phantom 
relative  to  the  plane  of  radiation  field  rotation  is  changed. 


Head  positioned  at  -20  degrees. 

Head  positioned  at  0  degrees.  Note  that  the  transverse  calculation 
plane  i3  coincident  with  the  plane  of  rotation  of  the  radiation 
field.  This  corresponds  to  the  plane  in  which  the  2-D  calculations 
were  performed. 


Head  positioned  at 
Head  positioned  at 
Head  positioned  at 
Head  positioned  at 
Head  positioned  at 
Head  positioned  at 
Head  positioned  at 
Head  positioned  at 
Head  positioned  at 
Head  positioned  at 
Head  positioned  at 


20  degrees. 
40  degrees. 
60  degrees. 
80  degrees. 
100  degrees. 
120  degrees. 
140  degrees. 
160  degrees. 
180  degrees. 
200  degrees. 
220  degrees. 


The  summation  of  dose  from  the  arcs  represented  in  these  views  is 
presented  as  the  3~D  rotation  dose  in  the  isocenter  transverse  plane  and 
the  isocenter  sagittal  plane  for  each  of  the  three  solar  flare  dose 
distributions. 
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